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We study the exclusive decays of B — > K^£ + £~ by the results in the perturbative QCD with 
the heavy quark effective theory and lattice calculations. We obtain the form factors for the 

0) 

transitions in the whole allowed region. Our predictions on the branching ratios 
of B -> 5 -> K*e + e-, and B -> K*fj,+fx~ are 0.53 ± 0.05ig;#, 1.68 ± 0.17±g;^, and 

U- 1.34 ±0.131 006 x 10 6 , where the errors are from the quark mixing elements and hadronic 

effects, respectively. We also find that our definitions of the T-odd observable and the up- 
down asymmetry of the K meson in B — > K*£ + £~ — > Ktt£ + £~ can be used to probe new 
physics. 



1 Introduction 



The flavor- changing neutral current (FCNC) processes of B — > K^H + £ (£ = e, //, r) are 
suppressed and induced by electroweak penguin and box diagrams in the standard model 
(SM) with branching ratios of e>(l(T 7 - 1(T 6 ) [1, 2, 3, 4]. Recently, the decay modes of 
B -> Kl+f (£ = e,/i) have been observed with Br(B -> i^+r) = (0.75±g;|f ±0.09) x lO" 6 
[5] and (OJSj^lotais) x 10 ~ 6 [ 6 > 7 ]> at the Belle detector in the KEKB e+e~ storage ring 
and the BABAR detector in the PEP-II B factory by using 29.1 fb' 1 and 77.8 fb' 1 data 
simples, respectively. Experimental searches at the B-factories for B — > K*£ + £~ are also 
close to the theoretical predicted ranges [5, 6, 7]. At BABAR, an excess of events over 
background with estimated significance of 2.8(7 has been observed and Br[B — > K*£ + £~) = 
(1.68±o:ii ± °- 28 ) x 10_6 has been obtained [7]. It is clear that these FCNC rare decays 
are important for not only testing the SM but probing new physics such as supersymmetric 
heavy particles in SUSY models, appearing virtually in the loop diagrams to interfere with 
those in the SM. 

It is known that one of the main theoretical uncertainties in studying exclusive hadron 
decays arises from the calculations of matrix elements. In our previous work [4] , we calculated 
the relevant transition form factors for B — > K(K*) decays in the perturbative QCD (PQCD) 
approach with the phenomenological wave functions, which are chosen by fitting with the 
experimental measurements of B — > Kir and B — > K*^ decays. Moreover, in order to 
compensate the increasing soft gluon effects in the slow recoil region, we used trial (in- 
dependent wave functions instead of unknown 6-dependent wave functions, where b is the 
conjugate variable of transverse momentum of valence quark. However, with the trial wave 
functions, it is difficult to estimate the errors in the B — > K(K*) form factors as well as the 
branching ratios and other physical observables in B — > K(K*)£ + £~ decays. 

To get a control of theoretical uncertainties in B — > K(K*)£ + £~ decays, in this paper, we 
first recalculate the corresponding form factors at the large recoil region of the momentum 
transfer q 2 — > with the improved PQCD approach that includes not only kr resummation, 
for removing end-point singularities, but also threshold resummation, for smearing the double 
logarithmic divergence arising from weak corrections [8]. The involved wave functions are 
used up to twist- 3, derived from QCD sum rule [9]. Then, we fit these results with those 
from the heavy quark effective theory (HQET) and the lattice QCD (LQCD) calculations in 
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the large q 2 region to get the whole q 2 allowed values. 

To emphasis the role of new physics, we define some P and T-odd observables. We 
will show that the effects of some observables are small in the SM but can be large and 
measurable in models with new physics. 

The paper is organized as follows. In Sec. II, we present the form factors of the B — > 
transitions in the framework of the PQCD. In Sec. Ill, we estimate the decay rates of 
B -> K^£ + £~ in the SM. We also study the polarization of K* in B -> K*£ + £~. In Sec. IV, 
we discuss P and T violating effects in B — > K*£ + £~ — > Krc£ + £~. We present our conclusions 
in Sec. V. 

2 form factors in B — ► K^*> transition 

To obtain the transition elements of B — > H (H = K, K*) with various weak vertices, we 
parametrize them in terms of the relevant form factors as follows:" 

(K(p 2 ,e)\V,\B(Pi)) = U(q 2 ){P,-^r-} + ^-fo(q 2 )q„ 

(K(p 2 ,e)\T^\B( Pl )) = M f ] {P-qq,-q 2 P,\, 

rriB + rriK L ' 

(K*(p 2 ,e)\V^\B( Pl )) = i V{ f ) e^ p e* a P^, 

m B + m K * 

(K*(p 2 ,e)\A^\B(p 1 )) = 2m K *A Q {q 2 ) e ^q^ + {m B + m K ,)A 1 {q 2 )[tl- e -^-q ll ) 

rriB + rriK* v q ' 

(K*(p 2 ,e)\T^\B( Pl )) = -miq^e^e^Py, 

(K*(p 2 ,e)\T*X\B(Pi)) = T 2 (q 2 )(e;P-q-e*-qP,)+T 3 (q 2 y-q(q,-^f-P,) (1) 

where V M = S7 M 6, A^ = S7^7 5 &, = sia^b, and T* v = sia^^b. To evaluate the q 2 - 
dependent form factors in Eq. (1), we use two QCD methods. One is for the large recoil 
region of small q 2 and the other for the zero recoil region of high q 2 . 

It is known that in the large energy transfer processes, the hadronic transition matrix 
elements can be calculated by the Lepage-Brodsky (LB) [10] formalism. However, the original 
LB formalism suffers logarithmic and linear singularities in twist-2 and twist-3 wave functions 
from the end-point region with a momentum fraction x — > 0, respectively. In order to 
"The relationships between various definitions of the form factors have been given in Appendix of Rcf. 

[4]- 
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handle these singularities, the strategy of introducing the kx resummation and threshold 
resummation has been proposed. It has been shown that due the induced Sudakov factors 
which make the \kr\ from A scale to (Am^) 1 / 2 with A = mg — mj [8] and m,b being the 
6-quark mass, the singularities do not exist in a self-consistent improved PQCD analysis 
[11]. Furthermore, due to the different properties of wave functions at end-point region, it 
is found that the power behaviors in B — ■> K(K*) form factors are the same in twist-2 and 
3 wave functions [12]. It shows the importance of twist-3 contributions. According to Ref. 
[12], we also know that other higher twist wave functions exist extra power suppression in 
n^o/ m B with m being the chiral symmetry breaking parameter so that they are neglected 
in our considerations. 

Recently, the applications of the PQCD approach to exclusive heavy B meson decays, 
such as 5 -> Kn [13], B -> nir(KK) [14, 15], B -> <pit(K) [16, 17], and B — > pK [18] decays, 
have been studied and found that all of them are consistent with the current experimental 
data. As known, in these two-body charmless decays, the squared momentum transfer is 
around m 2 M with M being a pseudoscalar or a vector meson. That is, the PQCD can work 
well at the region of q 2 ~ rn 2 M . Therefore, we will apply the predicted results of the PQCD 
to the large recoil region where the final outgoing meson M carries a large energy and 
momentum. As to the opposite region near q 2 \ max , we can use the relations among the form 
factor given in the HQET [19] by substituting the calculated values of the form factors from 
the LQCD [20] in them to get the remaining ones. Once the values of the form factors at 
both end edges of the allowed q 2 are determined, we can find the (^-dependent ones by fitting 
the forms 

FM 2 ) = ; - Fii °l 2 (2) 

1 + (TiS + (T 2 S Z 

where s = q 2 /m 2 B , Fj(0) are form factors at q 2 = 0, and a lj2 are the fitted parameters. 
We now show how to obtain the form factors at large and zero recoil: 

2.1 At large recoil 

In B — > Hi + i" (H = K, K*) decays, the B meson momentum pi, H meson momentum p 2 
and K* polarization vector e in the B meson rest frame and light-cone coordinate are taken 
to be 

Pi = ^=(1,1, 0±), P2 = ^^{V ,r H ,0 ± ), 
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*L = ^=^(^ 2 ,-r^,0 ± ), e T (±) = ^(0,0,l,±i) (3) 

where r\ ~ 1 — s and r# = m H /m B , while those for the spectators of B and /J sides are 
expressed as 

h = (o, x i~^f ' ^) » ^2 = (^2^|?7, 0, £ 2 ±^ , (4) 

respectively. In our calculations, we will neglect the small contributions from m u ^ s and A 
as well as m 2 H due to the on-shell condition of the valence-quark preserved. From the results 
of Ref. [9], the meson distribution amplitudes can be derived up to twist- 3 as follows: 

(tf(p)|5(^d(0),|0> = -^j\ xe ^{^ 5 Mx)+™^]Mx) 

(X*(p,e L )|s(z)^(0),|0) = £ dxe^imK^Mij^x) + Mk*{*) 

+m K *[I] lj (f) s K *(x)}, 

(lT(p, e T )|s(^(0),|0> = ^J^dxe^im^iM^^ + ^T Ml^x) 

+^^ie^ P ah5r}i^TP P n a ^K*(x)}, (5) 

P ' 71/ 

where n + = (1,0, Oj.) and n_ = (0,1, 0j_). In Eq. (5), (/>k(x), 4>k*{x) and <j>K*(x) are the 
twist-2 wave functions, and all the remaining ones belong to the twist-3. Their explicit 
expressions can be found in Ref. [9] . 

By using the LB formalism with including k F and threshold resummation, the form 
factors f + (q 2 ), f_(q 2 ), and fr(q 2 ) in B — > K can be written as 



U(q 2 ) = fi(q 2 ) + f2(g 2 ), 

W) = A(^ 2 )(i + 4)+/ 2 (g 2 )(i-4), (6) 



where 



pi roc r 

fi(q 2 ) = 8irC F m 2 B r K [dx] b 1 db 1 b 2 db 2 <Ab(zi, h) 4?k{x 2 ) ~ fticfa) 

Jo Jo 1 

xE(t^)h K (x 1: x 2 M:b 2 ), 
f2(q 2 ) = 8irC F m B / [dx] / b^b^d^ (j) B (xi, h) 

JO JO 

x{[(l + x 2 rj)(f) K (x 2 ) + 2r^((i -x 2 )<p t K {x 2 ) -x 2 (ff K (x 2 )) 
xE{t^)h K {x u x 2 MM) 

+2r K <P P K (x 2 )E(tM)h K (x 2 , x u b 2 , h) } , (7) 
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and 

pi roc 

f T (q 2 ) = 8nC F m 2 B (l + r K ) [dx] b 1 db 1 b 2 db 2 <l>B(xi,b 1 ) 

Jo Jo 

x{[Mx2) -r K x 2 {<t? K {x 2 ) -(t> t K {x 2 ))]E{t^)h K {x 1 ,x 2 ,b l M) 
+2r K( p p K (x 2 )E(t^)h K (x 2 , x u b 2 , h)}. (8) 

From Eq. (6), we find that /+(0) = /o(0). The evolution factor is given by 

E(t) = a s (t) exp(-S B (t) - S K (t)) , (9) 

where the Sudakov exponents S B (k) are given in Ref. [21]. The hard functions of h are 
written as 

h(x 1 ,x 2 ,b 1 ,b 2 ) = S t (x 2 )KQ( y /x^x 2 TimBbi) 

x[6(bi - b 2 )KQ( s fx^m B b 1 )lQ(^/x^firnBb 2 ) 

+9(b 2 - b^Koi^x^mBb^Io^y/x^qmBh)] (10) 

where the threshold resummation effect is described by [8] 
The hard scales t^ 1 ' 2 ^ are chosen to be 



= max(Jm B r)x 2 , 1/bi, l/b 2 ) 



t (2) = max(yjm%r)x 1 ,l/b 1 ,l/b 2 ). 

In Ref. [12], we have given detailed discussions and expressions for all form factors in the 
B — > K* transition based on the PQCD in the large recoil region. We emphasize that in the 
PQCD approach there is an identity at q 2 = 0, given by [12] 

A 2 (0) = (l + r K *) 2 A 1 (0)-2r K *(l + r K *)A (0), 



which is consistent with the leading order model-independent relation [1, 22, 23, 24, 25, 26] 

A 2 (0) = 11^^(0) - -^-A(0). 

1 — Tk* 1 — f K* 
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Table 1: Form factors for B -> K* in the LEET and PQCD with (I) u B = 0.40 and (II) 
uj b = 0.42. 





V(p) 


A (0) 


A 1 (0) 


A 2 (0) 


71(0) 


Ts(0) 


LEET [24] 


0.36 ±0.04 




0.27 ±0.03 




0.31 ±0.02 




PQCD (I) 


0.355 


0.407 


0.266 


0.202 


0.315 


0.207 


(II) 


0.332 


0.381 


0.248 


0.189 


0.294 


0.193 



In our numerical calculations, we use fs = 0.19 GeV, fx* = 0.21 GeV, = 0.17 GeV, 
m B = 5.28 GeV, m K * = 0.892 GeV, and c = 0.3 (0.4) for B -> X (X*), and we take the 5 
meson wave function as [21, 8] 



0_b(x, 6) = N B x 2 (l — x) 2 exp 
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2 , ,2 1,2" 






' 2 



(11) 



where o;b is the shape parameter [27] and Nb is determined by the normalization of the 
wave function, given by 



L 



fn 



(1x4>b(x, 0) = 



/o ^ v ' ' 2^2W C ' 

Since the shape parameter ojb and chiral symmetry breaking parameter m° K are free in the 
PQCD, in order to estimate the uncertainties, we choose (I) uo B = 0.40 and (II) uo B = 0.42 
for B -> K* and (I) uj b = 0.40 and m° K = 1.7 and (II) u B = 0.42 and m° K = 1.5 for 
B ^ K as the upper and lower bounds, respectively. From these values, we obtain the form 
factors of B -> AT at g 2 = as: (I) /+ = 0.354 and f T = 0.250 and (II) /+ = 0.303 and 
/t = 0.220, while those for the -B — > AT* ones are shown in Table 1. It is known that at the 
small q 2 region there exists a large energy effective theory (LEET) so that all form factors 
for heavy-to-light decays can be described by few independent functions [23]. In Table 1, we 
also show the results found by combining the LEET with the experimental data of B — > K*^ 
[24]. We note that our results for the form factors are different from those in Ref. [4], which 
shall be referred as (III). We remark that we have neglected nonlocal contributions from the 
four-quark operators. These effects can contribute in the 5 — 10% range to the form factors 
[28]. 
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Table 2: Form factors in B — > K from the HQET and lattice calculations at some large 
values of q 2 . 



q 2 GeV 2 


U(q 2 ) 


fo(q 2 ) 


f T (q 2 ) 


21 


1.61 ±0.11 


0.61 ±0.03 


1.54 ±0.12 


22 


2.05 ±0.12 


0.67 ±0.03 


1.99 ±0.13 



2.2 At zero recoil 

2.2.1 Form factors in B -> K 

According to the analysis of Ref. [19], under the heavy quark symmetry, the form factor of 
/t(<? 2 ) in B — > K can be written in terms of two independent form factors as 

Mg 2 ) = ^^(Utf) - f-(<i 2 )) (12) 

where f-(q 2 ) = (m B — m 2 K )(f (q 2 )— f+{q 2 ))/q 2 - Since there is no complete lattice calculations 
on the form factors in B — * K in the literature yet, the strategy to obtain /+ and /_ is that 
we utilize the relationships, which connect the form factors between B — > K and D — > K 
with the HQET, given by [19] 

f B( \ I ,m B 1/2 sa s (m b )^-(6/25)r m D m D D 

f + {V-p) = - )' {—, rl (1+ )f + (v P )-(l )f-(v-p) 

v, m.r, ^ rv _ ( m. . 1 / L m. n m.r> 



2 m B ^a s (m 

1 ,m B 1/2 / a s (m b )^-(6/25) 

2 m D ^a s (m c 



m B m B 



/£(„. P ) = L^'K^) 1(1 + -)/-<■« -—If (13) 



where ^a s (m;,)/« s (m c )j ^ ^ is the relevant renormalization factor. We note that the form 
factors f+(-)(v ■ p) are evaluated at the same value of v -p as f+,_Jv -p). In order to get the 
values at zero recoil, we adopt the lattice results in D — > K [29], where the predicted BR for 
D — > Kiv is consistent with the measurement [30]. From Eqs. (12) and (13), we find the 
form factors shown in Table 2 in the large q 2 region. In Table 3, we show the form factors 
in B — > K at q 2 = and fitted parameters of in Eq. (2). 

2.2.2 Form factors in B -> A"* 

We calculate the form factors in B — > A* by using the same technique as that in B — > K. 
Following the HQET, we have the relations given by [19, 31] 

2 _ m 2 B + q 2 -m 2 v V(q 2 ) m B + m v 2 

J iW J - ^ ; 1 o A n<7 J> 

2ms ms±my 2m B 
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Table 3: Form factors in B — > K at q 2 = and fitted parameters of a 1>2 in Eq. (2). 







/o(g 2 ) 


W) 


(I) q 2 = 


0.354 


0.354 


0.25 




-1.246 


-0.297 


-1.570 




0.251 


-0.400 


0.584 


(II) g 2 = 


0.303 


0.303 


0.220 




-1.229 


-1.212 


-1.496 


0"2 


0.219 


0.755 


0.492 



Tl(g) " T2(g) " -* - 2 2^ m s + m v " 2m B } 



mfj — m 2 / 



2m B mB + mv msq 

(m B + m v )A 1 (q 2 ) - (m B - m v )A 2 (q 2 )} . (14) 



m l + 9 2 - "^y 



2ms q 2 

We remark that the above identities are valid only for q 2 close to the zero recoil region where 
the HQET is reliable. We note that the relations in Eq. (14) are not complete as they mix 
terms of different order in l/m& in the zero recoil region [32]. The correct treatment has 
been given recently in Ref. [33]. From Eq. (14), we find that V(q 2 ), Ai(q 2 ), and T 3 (q 2 ) can 
be determined once T lj2 (g 2 ) and A 0t2 (q 2 ) are fixed. In Table 4, we display the form factors 
in B — > K*, where we have used Eq. (14) in the HQET and the lattice QCD results [20] of 
T 1}2 {q 2 ) and A 0j2 (q 2 ), which have been demonstrated to be consistent with the measurement 
in B — > K*^y. We remark that A (q 2 ) = Ti(q 2 ) has been taken in the lattice calculations. 
The form factors in B — > K* at q 2 = and fitted parameters of <J\ )2 in Eq. (2) are shown in 
Table 5. 

Table 4: Form factors for B — > K* from the lattice results of T 12 (g 2 ) and A 0t2 (q 2 ) with the 
HQET at some values of higher q 2 . 



q 2 GeV 2 


V(q 2 ) 


Ai(q 2 ) 


A 2 (q 2 ) 


Ti(q 2 ) 


Uq 2 ) 


Ts(q 2 ) 


16 


1.33 ±0.05 


0.44 ±0.02 


0.67 ±0.03 


1.14 ±0.04 


0.47 ±0.02 


0.67 ±0.03 


17 


1.50 ±0.05 


0.46 ±0.01 


0.72 ±0.02 


1.28 ±0.04 


0.48 ±0.01 


0.73 ±0.02 


19 


1.94 ±0.03 


0.49 ±0.01 


0.82 ±0.01 


1.66 ±0.02 


0.49 ±0.01 


0.86 ±0.01 
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Table 5: Form factors in B — > K* at q 2 = and fitted parameters of a lj2 in Eq. (2). 





W) 




A 1 (q 2 ) 




W) 


W) 


W) 


(I) q 2 = 


0.355 


0.407 


0.266 


0.202 


0.315 


0.315 


0.207 




-1.802 


-1.282 


-1.034 


-1.906 


-1.749 


-0.975 


-1.777 


0-2 


0.879 


0.249 


0.514 


1.168 


0.816 


0.632 


0.964 


(II) g 2 = 


0.332 


0.381 


0.248 


0.189 


0.294 


0.294 


0.193 




-1.721 


-1.228 


-0.829 


-1.801 


-1.671 


-0.721 


-1.677 


0-2 


0.744 


0.148 


0.166 


0.993 


0.684 


0.202 


0.794 



3 Differential Decay Rates and polarizations 

From the definitions of form factors in Eq. (1), in the SM the transition amplitudes for 
B — > K^£ + £~ (£ = e,fx) can be written as 

G F&\ t ([„eff, x , , 2 n , n s fr(q 2 ) 



M K 
and 

with 



2y/2n 



{[C e 9 ff (rif + (q 2 ) + 2m b C 7 ((i) 



Tn B + m K 



P^h^ + C 10 f + (q 2 )P,hV} (15) 



M 



(A) 



itie^e*»{\)P a qP + 6e*(A) + £ 3 e* • qP„ 



(16) 



(17) 



where we have set = (£ = e, fi) , X t = V tb V* ~ 0.041 ± 0.002 [34], Cf //} are the Wilson 
coefficients (WCs) and their expressions can be found in Refs. [4, 35], and i — 1 (2) for 
= (9j) w hh j = 1, 2, 3, defined in Appendix. 

For the decays of B — > K£ + £~, since the K meson is a pseudo-scalar, the only interesting 
physical observables are the decay rates themselves, whereas other observables such as the 
lepton polarizations, discussed in detail in Refs. [1, 2, 3, 4], are hard to be measured by 
experiments. The differential decay rates for B — > K£ + £~ are given by [4] 



dT K (s) G 2 F a 2 \X t \ 2 m 5 



B, 1 _ s) 3/2 



C e 9 ff (n)f + (q 2 ) + 2m b C 7 (n)f T (q 2 ) 2 + C w f + (q 2 ) *1 (18) 



ds 3 x 2 9 tt 5 

where s = q 2 /Mg. For B — > K*£ + £~ , however, one can study not only the decay rates and 
lepton polarizations but also physical observables related to the K* polarization, including 
longitudinal and transverse polarizations of K* and P and T-odd observables. To analyze 
the polarization of K*, we have to consider the decay chain B — > K*£ + £~ — > Ktc£ + £~, and 
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we choose the K* helicities as e(0) = {\p K *\,0,0, E K *)/M K , and e(±) = (0, 1, ±i, 0)/-\/2, 
the positron lepton momentum p t + = y/q 2 ^, sin 6>/, 0, cos 6i)/2 with E K * = (M| — M K * — 



g 2 )/2y / g 7 and = yi?^* — M^, in the q 2 rest frame, and the K momentum p K = 

(1, sin 9k cos0, sin^ sin 0, cos9k)Mk* /2 in the X* rest frame where denotes the relative 
angle of decaying plane between Kn and l + l~ . From Eq. (16), the differential decay rates 
of B — > K*i + i~ — > Kni + i~ as functions of angles 9 K and are found to be 

d cos 9 K d<j)dq 2 2 1 % 6 m 2 ? 1 J 

x{l6cos 2 ^ E l^°! 2 

i=l,2 



with 



+8sin 2 ^E (|^ + | 2 + |A4-| 2 ) 

i=l,2 

-8 sin 2 ^[ cos 20 E ReMtM~* + sin 20 £ ImMt M~* 

i=l,2 i=l,2 

+3tt sin2^[cos0( J ReAl?(A^+* -A^ 2 *) + i?e(A^+ 

+ sin0(/m7W°(7W+* + M 2 -*) -Jm(M+ + Mr)Al2*)]} (19) 



A<± = V9 2 ( ±2 b^lV9^i + 6), (20) 

where a = 1(2). We note that other discussions for various polarizations can be in Refs. 
[36, 37, 38]. From Eq. (19), by integrating all angles we obtain 



dT K *(s) G 2 F a 2 \\ t \ 2 \p\ 
ds 3 x 2 8 tt 5 



E EIMT, (21) 

A=+,0,-i=l,2 



where we have used that Br(K* — > Kir) = 1 [30]. 

The components M° a and Ai^ in Eq. (20) clearly denote the longitudinal and transverse 
polarizations, which can be extracted by integrating out the angle dependence in Eq. (19), 
respectively, and explicitly we have that 

+ sm 2 9 K E (\Mt\ 2 + \M;\ 2 )}, (22) 



=1,2 
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From Eqs. (21) and (22), we may define 

v ^ = roSwp - (23) 

, T(a = if^f^g , (24) 
as the normalized longitudinal and transverse parts with their ratio being 

? w ) = ^ / 2 n = ^ • ( 25 ) 



In Figures 1 and 2, we show the differential decay rates of B — > as functions of 

s = q 2 /Mg with and without resonant cc states, respectively. From both figures, it is obvious 
to see that the dilepton invariant distributions for the two modes are quite different. However, 
we expect that the distribution for B — > K£ + £~ , which contains only the longitudinal part 
due to the angular momentum conservation, is the similar to the corresponding part in B — > 
K*i + i~ . For B — > K*£ + £~, the transverse part, associated with 1/g 2 from the 7-penguin 
diagram described by CV(/i), gives the dominant contribution when q 2 goes to the allowed 
minimal values of Am 2 . This not only explains why the differential decay rates increase as 
q 2 — > but also indicates the reason for the difference between BRs in B — > A'*e + e _ and 

Using Eq. (21) and the fitted form factors in Tables 5 and 3, we obtain 

Br(B -> i^+r ) = (0.53 ± 0.05l°;J?) x 1CT 6 , 
Br(B -> K*e + e~) = (1.68 ± 0.17^) x 1CT 6 , 

5r(5 -> K*n + n~) = (1.34 ± 0.13+^) x 10~ 6 , (26) 

where the first and second errors are from A t and the hadronic effects, shown in Tables 3 
and 5, respectively. We remark that the central value of Eq. (26) for B — > K*e + e~ is the 
same as the BABAR recent measured value [7] and it is much smaller than that in (III) of 
Ref. [4] but that for B — > K*/j, + /j,~ a little larger. 

In Figure 3, we present the effects of K* longitudinal and transverse polarizations based 
on Eqs. (23) and (24). In Figure 4, we display the ratio of £(g 2 ) and we see that when 
s < 0.016 , where B — > K*^y is the main effect, and s > 0.339, where the longitudinal 
contribution is suppressed, the transverse part of the K* polarization becomes dominant. 
We find that both Vl and Vt as well as £(q 2 ) are insensitive to the hadronic uncertainties 
in (I) and (II). Moreover, the differences between the results of (I, II) and (III) are small. 
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4 T and P violating effects 

The terms with imaginary parts in Eq. (19) are related to T odd effects which, without 
final state interactions, are T violating and thus CP violating due to the CPT theorem. In a 
three-body decay, the triple correlations such as Sk-piX Pj are examples of the effects, where 
Sk denotes the spin vector carried by one of involving particles and p it j are the momentum 
vectors of outgoing particles. In the decays of B — > K*£ + £~, the spin Sk can be either the 
polarized lepton, S£, or the K* meson, e*(A). However, since the lepton polarization is always 
associated with the lepton mass and expected to be suppressed and less than 1% for the e 
or pi mode. On the other hand, the T odd effects with K* are free of the mass suppression 
and they can be large in models with new physics [12, 39]. To study these effects, we define 
[12, 39] 



r dT 

(Oi) = J O^u^ue^)— (27) 

where Ui(ug K ,ug (+ ) = ug K uo e+ /\u dK ug + | are sign functions with ug i being cos 6^ or sin#j. In 
the K* rest frame, we use the T odd momentum correlations as the operators in Eq. (27), 
given by 

n i- i (Pb ■ Pi+ x p K ){PB x p K ) ■ (p t + x p B ) 1 . 

°Tl = \Pb\ ^ v - 121- Z - 12 =oS m2 0, 28 

\Pb x PkV\Pi+ x pb\ 1 2 

rn i - i Pk-(pb* Pi+) . , {on , 

Ot 2 = \Pb\tz ^77^ —r = sm$, (29) 

\Pb x Pk\\Pb x pi+\ 

accompanied with sign functions of (sin 9k, sin 9 e+) and ut 2 (cos 9k, sin 9 e+), respectively. 
By defining the physical observables as 

* = M' m 

from Eqs. (27), (28) and (29) we obtain that 



E t =i,2 Im(MjM; 
4EA=+,o,-£ i= i,2|A^ 



„ Z^i=l,2- t " t V vl i -> vl i J / Q1 \ 

a ti = -— ^ nrix^ ( S1 ) 



3tt ImM^Mt* + MT) - Im(Mt + M^)Mf 

At2 ~ i6 e^Te^mF (32) 

We note that as shown in Refs. [12, 39] Atj. and At 2 depend on lmC% (/i)CV(//)* and 
ImCr{fJ>)Ci Q , respectively. In the SM, At 2 is zero since there are no absorptive parts expected 
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from C 7 and C w , whereas At x is non-zero due to that from Cg . However, it is clear that 
the T-odd observable of At 2 can be large in models with new CP violating phases beyond 
the SM. To illustrate a new physics result, in Figure 5 we show the T violating effect of 
At 2 as a function of s by taking two cases with the imaginary parts of WCs as follows: (i) 
ImC 7 (fi) = 0.25 and (ii) ImCj(fi) = 0.25 and ImCio — —2.0. One possible origin of having 
these imaginary parts is from SUSY [39] where there are many CP violating sources. Here, 
we only show the cases with (I) and (III) while those in (II) are almost the same as (I). It 
is interesting to note that the CP violating effect can be as large as 15% in both cases as 
shown in Figure 5. We remark that At x is much smaller than At 2 i n most of cases with new 
physics. 

From Eq. (19), we can also study another interesting physical observable associated with 
the angular distribution of sin 29 K cos in B — > K*£ + l~ — > Kirl + l~. This observable is 
defined as an up-down asymmetry (UDA) of the K meson due to its dependence of cos 6^. 
Explicitly, we define that 

n (£b x p K ) ■ (p e+ x p B ) 

Oud = —rz ^77^ =tt" = cos , 

\Pb x p K \\Pe+ x p B \ 

uj(ue K ,u de ) = u(cos0 K ,sm0e) , (33) 

as the operator corresponding to the UDA of K in the K* (— > Kir) rest frame and from Eqs. 
(19), (27) and (30) we find that 



a k < 2^ (Pud) _ 3n ReM^Mp - + Re(Mj - Mj)M° 2 * 

AuD{q } " dr/df ~ 16 E^Te^mF ' ( } 

which clearly violates parity. 

In Figure 6a, we show the Af} D (s) as a function of s = q 2 /M B based on the form factors 
given by the PQCD (I), (II) and (III) in the SM. It is interesting to point out that, as shown 
in the figure, in the SM A§ D (s) crosses zero point at s — 0.08 which satisfies the identity 

|2 



Re(h l9 * + h 2 g*) = -2 y ™ B *°J K * 1 Re(h 3 gl + h l9 *) . (35) 

We note that the point s is insensitive the QCD models but the Wilson coefficients of C 7 (/j,) 
and Cg(fi), especially the relative sign between them. To show the result, in Figure 6b we 
present two extreme cases of (i) Cr(fJ>) = —C-j{ii)sm (solid curve) and (ii) Cg(/x) = —Cg(/j,)sM 
(dash-dotted curve) with the remaining Wilson coefficients the same as those in the SM, 
respectively. From the figure, we see that the distributions in the two cases are quite different 
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and moreover, s disappears, i.e., it exists only if CV(/x) and C 9 (/x) have an opposite sign 
since Cr(fj)sM < and Cg(/i)sM > 0. It is clear that s provides us a good candidate to 
explore new physics due to the insensitivity of the QCD models and dependence on the WCs. 



5 Conclusions 

We have studied the exclusive decays of B — > K^£ + £~ by the results in the PQCD with 
the HQET and LQCD. We have given the form factors for the B — > ifW transitions in the 
whole allowed region, which are consistent with those from other QCD models. We have 
found that the branching ratios of B — > K£ + £~, B — > K*e + e~, and B — > K*fi + fi~ are 
0.53 ± 0.05l8;$, 1.68 ± and 1.34 ± 0.13lg;JJ x 10~ 6 , respectively. We have shown 

that, in B — > K*l + l~ — > Krc£ + £~, the T-odd observable of ^2 which is unmeasurably small 
in the SM could be as large as 20% in models with new physics, while the P-odd up-down 
asymmetry of A§ D (s) vanishes at sq ~ 0.08 in the SM but it behaves quite differently if new 
physics exists, which provide us unique probes of non-standard physics. 
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Appendix 



The parameters of ^1,2,3 an d (71,2,3 in Eq. (17) are defined by 

h = cm v{q2) + 2 -^(n)CMT 1 ( q % 



ttib + m K * q 



h 2 = -CMirnB + m^A^-^P-qCMniq 2 ), 



h = CM A f ] +^ 7 (^)(T 2 (g 2 ) + -^T 3 (g 2 )), 
mR + m^» a z v P ■ a 1 

9i = C w - 



tub + rax* q z F-q 
V(q 2 ) 



m B + m K * 



92 = -C 10 (m B + rriK^A^q ) , 
fl's — ^10 ; • 

TUB + VtlK* 
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Figure 1: Differential decay BR of B — > Kfi + fi~ with (a) and (b) representing the results 
with and without resonant effect, respectively. The solid (dash-dotted) curve stands for the 
upper (lower) bound with the allowed region being shaded. 




Figure 2: Same as Figure 1 but B — > K*fi + fi + . 
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Figure 3: Normalized Longitudinal (a) and transverse (b) polarizations in B — ► K*fi + fi~. 
The solid, dash-dotted and dashed curves denote the results from the PQCD (I), (II) and 
(III) respectively. 
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Figure 5: T violating effect of At 2 (s) for (i) imCV(/i) = 0.25 (solid curves) and (ii) 
lmCi{\x) = 0.25 and ImCi = —2.0 (dash-dotted curves), where the bold and thin lines 
correspond to the PQCD (I) and (III), respectively. 
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Figure 6: A§ D (s) as a function of s. In (a) the solid, dash-dotted and dashed curves represent 
the SM contributions based on the form factors in the PQCD (I), (II) and (III), while in (b) 
the solid (dash-dotted) and upper (lower) dashed curves are for CV(/z) > (Cg(fi) < 0) in 
(I) and (HI), respectively. 
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